Instability of planar vortices in two-dimensional easy-plane Heisenberg model with 

distance-dependent interactions 
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It is known that magnetic vortices in two dimensional Heisenberg models with easy-plane 
anisotropy exhibit an instability depending on the anisotropy strength. In this paper, we study the 
statistic behavior of the two-dimensional easy-plane Heisenberg models with distance-dependent 
interactions, J xy (r) and J z (r) for in-plane and out-of-plane components. We develop analytical 
and numerical methods for accurate determination of critical anisotropy, above which out-of-plane 
vortices are stable. In particular, we explore the vortex formation of the Gaussian-type interaction 
model and determine the critical anisotropy accurately for square, hexagonal and triangular lattices. 



I. INTRODUCTION 

The magnetism in two-dimension has been the subject 
of continued interest in the last few decades, which is 
aroused further due to recent additions of new materials 
such as Cu-based high-temperature superconductors and 
advances in both numerical and experimental capabili- 
ties. The classical two-dimensional anisotropic Heisen- 
berg model (CTDAHM) with the spin Hamiltonian, 



H 



-K 



( Q x Q x 



(1.1) 



for a coupling constant K, has been one of the most 
studied models for magnetic systems. For A = 0, it be- 
comes the so-called XY model, which is well understood 
following the work of Berezinskiii and Kosterlitz and 
Thoulcss 2 . It has a phase transition at Tkt, named the 
Kosterlitz-Thouless phase transition, which can be char- 
acterized by a vortex- antivortex unbinding*. For A < 1, 
it belongs to the universality class of the XY model, and 
it has been known that there are two different types of 
nonlinear solutions, called "in-plane" and "out-of-plane" 
vortices*. There have been several attempts with nu- 
merical simulations and analytic calculations to obtain 
the critical anisotropy A c , above which the out-of-plane 
vortices are stabilized. At first, Gouvea et air* inves- 
tigated the development of the out-of-plane vortex by 
using a combined Monte Carlo-molecular-dynamics tech- 
nique. They found that initial out-of-plane vortex relaxes 
to a planar one below the critical anisotropy A w 0.72, 
0.86, 0.62 for square, hexagonal and triangular lattices 
respectively. Costa and Costal obtained a more precise 
measurement A c ~ 0.709 for a square lattice by fitting 
the out-of-plane squared component (S z ) 2 at the cen- 
tral peak as a function of A. Wysin made an analysis 
of the core region of a vortex on a discrete lattice and 
obtained an estimate of A c « 0.7044 for the square lat- 
tice*, which is improved further to a very accurate value 
of A c « 0.703409 by iteratively setting each spin's xy 
components to point along the direction of the effective 
field due to its neighbors. 

Recently, we showed that the pattern formation in pri- 



mary visual cortex is also related with the vortex dy- 
namics in magnetism 8 . It is observed in vivo that the 
development of strong scalar components (ocular dom- 
inance column) near the singularity centers at orienta- 
tional components (orientation column). There are at 
least three different types in observed ocular dominance 
patterns in vivo, which correspond to the Ising type, sta- 
ble out-of-plane vortex, and stable in-plane vortex be- 
havioral. The determination of the critical anisotropy is 
important in describng the region of different behaviors 
in vortex patterns. 

In this paper, we consider the two-dimensional easy- 
plane (XY) symmetry Heisenberg model with distance- 
dependent interactions : 

H = - {J*v(n 3 )(S*S* + SfS*) + J z ( rij )SfS]}(1.2) 

for the classical spin variables Si = (Sf , Sf , Sf) and the 
distance dependent exchange energy J M (r) = e M / A1 (r)/2 
(/i = xy or z). The dimcnsionless neighborhood interac- 
tion function I^ir) can have a various shape but is ex- 
pected to be a smooth function and approach zero as 
r — > oo. Comparing Ea. l|1.2|l with Ea. (|l.lfl . we may 
expect that the anisotropy in the strength of interac- 
tions between in-plane and out-of-plane components is 
A — Jz/Jxy and obviously A = s z /e xy if I xy {r) = I z (r) 
for all r. But a more clear definition of the anisotropy 
parameter A will be needed for the case of arbitrary in- 
teraction functions with I X y{f) ^ Iz(f). Takeno and 
Homma showed that the anisotropic Heisenberg model 
with distance-dependent interactions has two different 
minima in H depending on the ground-state energy of 
in-plane and out-of-plane components*^. We determine 
the anisotropy A by the ratio of the ground-state energy, 
which is the boundary of the crossover behavior between 
the XY and Ising models. The critical anisotropy now 
may depend on the shape and parameter values of Jy,(r). 

In Sec. m we build an effective Hamiltonian in a con- 
tinuum approximation for the arbitrary exchange energy 
function J(r). The continuum theory helps to predict 
or understand the major behavior of spin configurations 
and vortex formation. But the exact estimation of the 
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critical anisotropy cannot be obtained by a continuum 
approximation because of the singularity near the vortex 
corei Sometimes two distinct spin models can yield the 
same Hamiltonian in a continuum approximation. We 
will show that different values of the critical anisotropy A c 
is possible in spite of the common approximated Hamil- 
tonian (Sec. Ej. In Sec. EH we generalize and extend 
Wysin's methoc&ii, which is the analytic method for 
obtaining the exact critical anisotropy A c in the distant 
neighbor interaction models. We clarify our method by 
calculating critical anisotropies of the CTDAHM with 
7(r) = 6(r — 1) for different lattice types and compar- 
ing them to values. In Sec. IIVI we exhibit a simula- 
tions! method for the determination of the critical value 
of anisotropy above which out-of-plane components de- 
velop. We find that the out-of-plane components start 
to develop for A > A c where A c are 0.7035, 0.8330 and 
0.6129 in four-digit accuracy for square, hexagonal and 
triangular lattices, respectively. These values from sim- 
ulations are the closest values to the analytic calcula- 
tions^. In Sec. we investigate the critical anisotropy 
in the two-dimensional easy-plane Heisenberg model with 
Gaussian type interactions of i M (r) = exp(— r 2 /2cr 2 ) and 
the anisotropy A = (£ z /s xy )(cr z /o'l y )- We obtain the 
critical anisotropy A c for A = s z /e xy (cr xy = a z ) and 
A = u 2 /a xy {e xy = e z ) depending on the interaction 
ranges tJ M for different lattice types and find a general 
behavior of A c . 

II. EFFECTIVE HAMILTONIAN APPROACH 



fourier transform of the expansion in Eq. 12.. '{1 : 

H ~ -NJ S S 2 + ^- J 'd 2 rj a (l-m 2 )(V0) 2 

+a[l - 6' (I - m 2 )] ,[ Vm) * + 4£m 2 l (2.5) 
(1 — m z ) J 

where J s = J xy {0), J p = -j£ y (0)/a 2 , 5 = 1- 

j,(0)/J*y(°)> S ' = 1 - j"(0)/j"y(0), a = 2J p /J s and 
rn(r) = S z (r)/S for the nearest-neighbor distance a. The 
exchange energy J^(r) need not be fully positive, how- 
ever, in ferromagnetic behavior systems. Comparing with 
the approximated Hamiltonian of the CTDAHM in con- 
tinuum theory^, we expect a ferromagnetic behavior if J s 
and J p are positive. The anisotropy A is described as 

A= 1-5 = J Z (0)/J xy (0), (2.6) 

which is the ratio of the ground-state energy for in-plane 
and out-of-plane components. Eci. (|2.5|l has two differ- 
ent solutions depending on the anisotropy in A. The 
total ground-state energy in Ea. (|2.5|) . (H)o, is given 
by -NJ xy {0)S 2 (or -NJ z (0)S 2 ) and m 2 (r) = (or 
m 2 (r) = 1) for all r when A < 1 (or A > 1), which 
is the XY (or Ising)-like solution as shown by Takeno et 
alM* The anisotropy satisfies A = e z /e xy if I xy (r) = I z (r) 
for all r. However ^ Iz{f), the anisotropy can be 

expressed as A cx (<r z /a xy ) d for the interaction range 
in the d-dimensional lattice because the exchange energy 
is proportional to the number of interaction pairs in fer- 
romagnetic systems. 



The classical spin vector, Si = (Sf,Sf,Sf), can be 
specified by two angles of rotation </>j and 9i, 

Si = S(sm9i cos <pi, sin 9i sin <^>j, cos 0i). (2-1) 

In the momentum space representation, we write 

i 

where the vectors q are restricted to the first Brillouin 
zone of the simple cubic <i-dimcnsional lattice. Substitut- 
ing Eas. (|2.1|l and l|2.2|) into the Hamiltonian in Ea. l|1.2|) . 

we get 

h = -J2{ J M(s q J^+s qy s^ qy ) 

q 

+J z (q)S qz S- qz } (2.3) 

where 

J M (q) = ^J M (r) e -^ r (2.4) 

r 

for fi — xy or z. If J^q) has a maximum point at q = 0, 
the approximated Hamiltonian is obtained by the inverse 



III. ANALYTIC CALCULATIONS OF 
CRITICAL ANISOTROPY 

The Hamiltonian in Ea. (|1.2|l . can be rewritten as 

H = -^2{Jxy(Tij)\Jl - m?jT - TO 2 COS(0 4 - 4>j) 

+J z (r ij )m,im j }, (3.1) 

where 4>% = tan _1 (5,f / Sf ) and m, = S* /S. For arbitrary 
\rrii\ < 1, we get 

j 

-J z {r ij )m j } = Q. (3.2) 

This always has an in-plane vortex solution, to^ = for 
all sites i. For the out-of-plane solutions, which have 
non-zero rrii components, the determinant of the matrix 
W should vanish, where 

W a f} = S aj3 ^ E J xy {Tij)cOs{4>i - <j)j) 

-EE J *( r «)- ( 3 - 3 ) 
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FIG. 1: The vortex center is located at (0,0) and the sites 
are presented by different symbols according to the distance 
from the vortex center, (solid squares for r = l/y2, solid 
triangles for r = y/W/2, etc.) 



For an arbitrary exchange energy function J^(r), /i = xy 
or z, the critical value of the anisotropy can be deter- 
mined by a root finding method for det|W(A)| = 0. 
When I xy (r) — Iz(r) for all r, so A = e z /e xy , the bound- 
ary for the existence of the out-of-plane solutions can be 
obtained from the eigenvalue problem, 



for 



and 



A a m a — Xwapmp = 



(w a p/A a )mi3 = (l/A)m Q 



E E 7 ^) 00 ^ 



W a (3 



i M I ^ E J ( r «) 



E 



(3.4) 



(3.5) 



(3.6) 



(3.7) 



where a or /? is a block index rather than a site index 
when there exists a symmetry in the lattice geometry and 
M a = { i | ri — r a } are sets ordered by the distance from 
the vortex core (1 < a < L). The maximum eigenvalue, 
which is equal to 1/A C , is determined by using only the 
nearest components to the core center, mi,rri2, ■ ■ .mjj 
(1 < Af < L, where Af is the number of the nearest 
components). 

Note that we assume that the center of the vortex is 
located at the origin of a coordinate system. The in-plane 
angles are chosen to be the usual solution of the in-plane 
vortex in the continuum theory, <fii — tan -1 (yt/xi), or 




FIG. 2: The critical anisotropy A c as a function of \og 2 Af 
for the anisotropic Heisenberg model with Gaussian type in- 
teractions for a = 2,4 and 8. As a increases, the size of the 
out-of-plane vortex core increases. 



For the CTDAHM in Ea. lfl~T|) . the interaction function 
is given by I(r) — S(r — 1) and 



At 

A 2 
w n 

U>21 



4 
4 

7E 

W22 



(1 



v/5 VT3' 



(3.9) 



where A'h = {±(1/2, 1/2), ±(1/2, -1/2)} and M 2 = 
{±(3/2, l/2),±(3/2, -l/2),±(l/2, 3/2),±(l/2, -3/2)} 
for the square lattice in Fig^ In this case, w a p is just 
the number of the nearest neighbors that belong to the 
set Mfj for rrij € M a . The critical anisotropy, where the 
2x2 determinant vanishes, is given by 



/V 



A X A 2 



A 1 + 2A 2 



0.7157 



(3.10) 



/) 



(3.8) 



which agrees well with the results of Wysini. Wysin also 
showed the A c is determined by the zero of a 3 x 3 de- 
terminant as 0.7044. We obtain a smaller critical value 
of A c w 0.6974 for Af = 4, which converges to 0.6941 
in four-digit accuracy for Af > 16. This critical value 
show a deviation from the result in numerical simula- 
tions. Wysin showed that there are some differences in 
the discrete solution from the continuum results espe- 
cially near the vortex coreii. Similar to Wysin's process, 
we eliminate this difference between discrete and contin- 
uum solutions by evolution of each spin's xy components 
to point along the direction of the effective field due to 
its neighbors, and we obtain A c « 0.703420, which is very 
close to the value obtained by Wysin's, A c 0.703409, 
and one from our simulations in Sec.|^(A c ra 0.7035). 

For the nearest-neighbor interactions in the CT- 
DAHM, the size of the vortex core is small, which leads 
to precise determination of the critical anisotropy for rel- 
atively small Af. But as the interaction range increases, 
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the size of the vortex core increases and more compo- 
nents of the matrix are needed to obtain the accurate 
value of A c . Fig. [3 shows the computed anisotropy A c 
as a function of Af for different interaction ranges c M in 
Gaussian type interactions. The convergence of different 
lines for different ranges gives an estimate of the vortex 
core size r v . For the average number of the degeneracy g 
per block M a , gj\f c means the number of spin sites which 
contributes to the determination of A c . The radius of the 
vortex core can be estimated as 



(3.11) 



with g approaching 8 for large j\f c for the square lattice 
or the D4 symmetry geometry. 




(a) 



IV. SINGLE- VORTEX SIMULATIONS 

In general, it is quite difficult to determine the criti- 
cal anisotropy precisely with a simulated annealing ap- 
proach using Monte Carlo simulations. In particular, for 
the distance-dependent interaction models, the longer an- 
nealing time is required due to the larger lattice size and 
the larger cutoff range. Our numerical simulations are 
carried out by evolving dynamical equations as a gradi- 
ent flow, which allows a more accurate determination of 
the vortex behavior near A ~ A c . The discrete equations 
of motion used in simulations are 



dt 



dt 



OH 

d(j>i 

-2 J xy (nj ) sin 9i sin 9j sin( 

i 

dH 
"d~9~. 



<t>j) (4- la) 



= 2 ^ { Jxy {nj ) cos Oi sin 9j cos(<fo - <j) ) 
3 

— Jz ( r ij ) sin Oi cos Oj } (4. lb) 

at zero temperature. The simulations start from an ini- 
tial state of a planar vortex on the lattice center and 
small random fluctuations in out-of-plane components 
(|5,|<10- 6 ). 

A problem in the single vortex simulations is the effect 
of the area boundaries. The distribution of spin compo- 
nents is effected by the area boundary as shown in Fig.|3j 
At first, the contour lines of iso-phase are stretch out 
radially from the vortex core. But after dynamical evo- 
lution, they tend to meet perpendicularly with the area 
boundary due to the equilibrium condition SH/Scj) ~ or 
V 2 ~ 0. If the lattice size is not sufficiently large, the 
distribution of spin components near the vortex core is 
effected by the area boundary. For the spin system with 
distance-dependent interactions, the size of the vortex 
core increases, so that the larger lattice size is necessary 
for simulations. Note also that the shape of the lattice 




(b) 




(c) 

FIG. 3: Spin components near in-plane and out-of-plane vor- 
tex formed in the easy-plane Heisenberg model with Gaussian- 
type interactions, (a) The Initial state with an in-plane vor- 
tex at the lattice center with small random fluctuations of 
\S Z \ < 10 -6 and the states after evolution for (b) A = 0.72 and 
(c) A = 0.74 (A = e z /e X y, Oxy — <r z — 1, square lattice type 
and 50 x 50 lattice size). In-plane spin components are rep- 
resented by contour lines where <f>i = nir/8 for n = 0, 1 ... 15. 
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is matter more than its size. Note that the single vortex 
is attracted to the area boundary like an electric charge 
near the conductor surface. If there is an asymmetry in 
the lattice, the centered vortex may move along a cer- 
tain direction. An asymmetry in the shape of the cutoff 
range can also induce the wandering of the vortex core. 
Especially the single-vortex simulations for the triangu- 
lar lattice turned out to be more difficult and sensitive 
to the shape of the lattice and the cutoff range. 

We measur e the growth of the out-of-plane components 
by yl ((S z ) 2 ) . The measurement of the squared compo- 
nent {S z ) 2 at the central peak can be another choice. 
However, tracing single or few components may fail to 
determine the vortex stability since the vortex core may 
start to move. In case that the vortex wanders, the fluc- 
tuations can be measured by \J ((S z ) 2 ). To clarify the 
simulation results, we confirm the monotonical increase 
or decrease after the inflection point. 

Our simulational method is very fast, yielding the com- 
puted critical anisotropy that is more accurate than those 
by previous approaches. For example, Fig. 0] shows the 
results of dynamical evolution with the CTDAHM for 
the square lattice. We find a monotonical growth (or 
decay) of out-of-plane components after some time for 
A = 0.7035 (or A = 0.7034, 0.7033). This graph sug- 
gests that the critical anisotropy lies between 0.7034 < 
A c < 0.7035, which agrees well with analytic calculations 
by Wysin in Ref. 9 (A c « 0.703409) and ours in Sec ITU 
(A c « 0.703420). Similarly, we obtain A c w 0.6129 for 
the triangular lattice and A c « 0.8330 for the hexagonal 
lattice, which also agree well with the results of Wysin 
(A c « 0.612856 and 0.832956 respectively^) . 




(c) 

FIG. 5: Deviation in the discrete solution from the continuum 
result near the vortex core for the square lattice, where the al- 
and y-axis are the lattice coordinates and the z-axis denotes 
)| for (a) the nearest-neighbor interactions 
(b) the Gaussian type interactions with a xy = 1 and (c) the 
Gaussian type interactions with a xy — 4. 

V. EXPERIMENTS IN GAUSSIAN TYPE 
INTERACTIONS 

As a typical example of the spin systems with an arbi- 
trary distance-dependent interaction types, we study the 
model with Gaussian type interactions where 

/^r) = exp(-r 2 /2<7 2 ). (5.1) 

In the continuum limit, the fourier transformed interac- 
tion, 

/„(«)=7r(^) a exp(-ay/2), (5.2) 

has a maximum at q = and the approximated Hamil- 
tonian can be written as Ea. l|2.5|l with J s = ne xy a xy / a 2 , 
5=1- {e z /e xy )(a 2 /a 2 xy ), §' = 1 - (e z /e xy )(af/a xy ), 
a 2 = 2a 2 and the anisotropy 
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FIG. 6: The critical anisotropy A c as a function of a xy in the 
anisotropic Heisenberg model with Gaussian type interactions 
for (a) A = Ez/exy {oxy = a z) and (b) A = ff*/o^ B (E xy = e z ). 
The critical anisotropy A c are determined by the determinant 
of the 1024 x 1024 matrix. The symbols denote the results of 
the numerical simulations. 

The critical anisotropy A c is not uniform but varies de- 
pending on e M and . We can investigate the behaviour 
of the critical anisotropy from Gaussian type interactions 
for two different general cases; the exchange strengths are 
proportional to each other, J z (r) / J xy (r) = const for all r 
or A = s z /e xy , and the exchange energy are proportional 
to a 2 L for the interaction range in the two-dimensional 
system or A = o%/o% y . 

In analytic calculations, a large number of the near- 
est components, TV, is needed for the convergence of the 
critical value A c as shown in Fig. [21 The critical values 
of the anisotropy are determined for Af — 1024. For the 
case of A = e z /e xy , the critical anisotropy is found by the 
condition that the determinant of the J\f xj\f matrix, W, 
vanishes where 

w a0 = 5 a p E^p^^/ 2 ^) 00 ^-^) 
- x E E ex p(-4/ 2 0- ( 5 - 4 ) 

The maximal eigenvalue A of Ea. H5.4fl gives the critical 
anisotropy A c as in Eq. (|3.5f) . The angles of in-plane com- 



TABLE I: Linear fitting of a xy \/\: as a function of a x 



Lattice 
type 


A — E z 1 Exy 




slope y-intercept 


slope y-intercept 


Square 
Hexagonal 
Triangular 


0.764928 0.082300 
0.765017 0.116011 
0.764967 0.068548 


0.625910 0.181816 
0.629297 0.243578 
0.624828 0.154744 



ponents 4>i are initially set to their continuum vortex val- 
ues, 4>i = tan -1 (yi/xi), for each site i and evolved until 
they approach the equilibrium state according to the dy- 
namical equation, 




= -2 V J xy {rij)sin{<j)i - fc) (5.5) 



j 

from Ea. H4.la)) with m t = for all sites, that is, the dis- 
crete dynamics in the absence of out-of-components or 
the pure XY model limit. Note that the spin configu- 
ration of the system with distance interaction, however, 
corresponds well to the continuum results near the vortex 
core in contrast to one with the nearest-neighbor interac- 
tion. There are much differences between our results and 
the classical models with the nearest-neighbor interac- 
tions, which decreases as the interaction range increases 
for the Gaussian type interactions. Fig. [3] shows the 
differences between the discrete solutions and the con- 
tinuum solutions. The discrete solutions are obtained 
after the evolution from the continuum solution fa — 
tan (j/j/arj) by Ea. H5.5f) . avoiding the area boundary ef- 
fect by using a round shaped lattice boundary (R = 80). 
However, we solve the equations with the initial angles 
or setting cos(c6; — 4>j) = ?i • fj for Gaussian type in- 
teractions, and which leads to the accurate results. For 
A = a 2 /a 2 y , the matrix W is given by 

W aP = 5 af 3 E exp ^~ r ?)/ 2cr ^) cos (^ _ 

ieM a j 

-EE ex p(-4/ 2A <4)> ( 5 - 6 ) 

which is difficult to be converted as an eigenvalue prob- 
lem as before. The critical anisotropy A c is determined 
by /(A) = det |W(A)| in Ea. H5.tif) using the root finding- 
method, which is in agreement with the results of nu- 
merical simulations within the accuracy of two- or three- 
digits. We show the analytic and simulational results for 
square, hexagonal and triangular lattices in Fig. |SJ 

The obtained critical value of anisotropy decreases as 
the interactions range a decreases. Note that the com- 
puted values of the critical anisotropy A c are in the order 
of hexagonal, square and triangular lattices type for both 
A = e z /e xy and A = a 2 /a 2 y as in the CTDAHM. Fig.0 
shows the plot of o xy \f\, as a function of a xy , showing 
a linearity between them. The slopes for different lat- 
tice types can be grouped as two classes with A = e z /e xy 
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FIG. 7: OxyVK: as a function of a xy . The upper branch 
corresponds to A = e z /e xy and lower branch to A = a^/a^y. 
When A = a^/cr^, y , y-axis means a z . 



and A = a\ja\ v (Tabled}. If the linearity is maintained 
for all a xy , the asymptotic value of \f\^ will be converge 
to the slope of the line for large a xy . With the help of 
the asymptotic solutions in Ea. (|5.4l) and Eq. H5.6f) . we can 
guess that the critical anisotropy A c will converges to zero 
and the slope of lines in Fig. will become smaller for 
large a xy 

When e xy = e z and a xy = a z — l/v2> the approx- 
imated Hamiltonian with Gaussian type interactions in 
Eo. (|2.5|) becomes that 



H ~ -NJ s S 2 + ^f- J rf 2 r|(l-m 2 )(V0) 2 

+ [l-,5(l-m 2 )] 7 ^^+4 ( 5m 2 l, (5.7) 

(1 — 777 J J 



which is equivalent to that of the CTDAHM 5 with a cou- 
pling constant of nearest-neighbor interactions J s /2. In 
this case, they have the common ground-state energy and 
the common excitation energy due to formation of a vor- 
tex for the same lattice type and lattice size. However, we 
find that the critical anisotropy A c in Gaussian type in- 
teractions is 0.7938 from analytic calculations and 0.7941 
from simulational experiments for square lattice, which 
are different with A c w 0.7035 in the CTDAHM. These 
results mean that the solutions of the critical anisotropy 
from the continuum limit approach are inadequate. 
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FIG. 8: Crossover behavior for interaction ranges. 



VI. DISCUSSION 

We have studied the vortex formation in the two- 
dimensional anisotropic Heisenberg model with distance- 
dependent interactions and determined the critical 
anisotropy A c precisely for Gaussian type interactions. 
The continuum theory helps to predict how spin config- 
urations will develop for an arbitrary shaped exchange 
function J M (r). If the exchange energy function Jn(r) 
satisfies certain conditions, the approximated Hamilto- 
nian in the continuum limit hows a similarity to that of 
the CTDAHM, which leads to the similarity in two types 
of static vortex solutions. However, the continuum the- 
ory breaks down in determining the critical anisotropy 
as noted by WysiiA We also give an example that two 
different anisotropic Heisenberg models with a common 
approximated form in the continuum limit can have dif- 
ferent critical anisotropies. 

It is very difficult to determine the critical anisotropy 
for an arbitrary type of the exchange energy function 
Jjj,(r). In this paper, we have studied two typical cases by 
using Gaussian type interaction models. One case is that 
Jxy(r) and J z (r) are proportional for all r and another 
case is that J xy (r) and J z (i") have a common distribu- 
tion type with different interaction ranges. The results 
for these cases suggest some general behaviors of A c . The 
critical anisotropy is the smallest for the triangular lattice 
and is the largest for the hexagonal lattice among three 
lattice types explored. This is also the case for the clas- 
sical model because the spins near the vortex core tend 
to be more parallel, reducing their total exchange energy 
further, so that the lattice type with more nearest neigh- 
borhoods can reduce the out-of-plane exchange energy 
further while increasing the in-plane exchange energy for 
lower value of A. Fig. [H] shows the crossover behavior of 
spin dynamics when the anisotropy between the in-plane 
and the out-of-plane energies is manifested only by the 
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difference in interaction ranges. If the interaction range 
of the in-plane is larger than that of the out-of-plane 
(&xy >oi), the number of exchange pairs in the in-plane 
case is larger for the ferromagnetic system and the system 
belong to the same universality class as the XY model. 
The region where the out-of-plane vortex is stable can be 
approximately found in the region with a < a z /a xy < 1, 
for a reasonable range of a. If a xy is very large in com- 
parison with the nearest-neighbor distance, a nearly zero 
critical anisotropy may occur. 



Our study for the precise determination of the criti- 
cal anisotropy is motivated from the problem of the map 
formation in the cerebric cortex rather than magnetism 
systems. The results for the crossover among the regions 
with Ising, stable out-of-plane or stable in-plane vortex 
behaviors depending on the interaction ranges in Fig. |H1 
correspond well with the data from the animal exper- 
iments for different types of the map formation in the 
visual corte»i. 
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